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Summary 

This  paper  is  an  attempt  to  understand  and  measure  dispersion 
effects  of  factors  in  factorial  experiments.  The  simplest  setting  is 
considered  in  order  to  develop  better  comprehension  and  insight.  The 
properties  of  the  proposed  descriptive  measures  are  examined.  A  method 
of  "adjusting"  residuals  and  its  use  in  measuring  dispersion  effects  are 
discussed.  Illustrative  examples  are  also  given.  The  problem  considered 
in  this  paper  arises  in  quality  control  studies  and  the  methodologies  are 
applicable  to  industrial  experiments. 
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1 .  Introduction 

An  important  problem  in  quality  control  studies  is  to  find  an 
optimum  combination  of  levels  of  control  factors  in  achieving  stability 
against  noise  factors  (see  Taguchi  and  Wu  1985).  Both  "location"  and 
"dispersion"  effects  of  factors  are  pertinent  to  measure  from  the  data  in 
resolving  this  problem.  This  article  considers  the  problem  of  measuring 
dispersion  effects  of  factors  in  both  replicated  and  unreplicated 
factorial  experiments.  The  concept  of  dispersion  effects  in  factorial 
experiments  was  considered  in  the  work  of  G.  Taguchi  (see  Taguchi  and  Wu 
1985)  for  replicated  factorial  experiments  and  in  the  work  of  G.E.P.  Box 
(see  Box  and  Meyer  (1986))  for  unreplicated  factorial  experiments. 
Factorial  experiments  may  be  complete  or  fractional  factorial  under 
completely  randomized  designs.  Although  for  clarity  we  consider  2ra 
factorial  experiments  in  this  article,  the  ideas  presented  can  easily  be 
generalized  to  any  symmetric  or  asymmetric  factorial  experiments.  Kackar 
(1985),  Phadke  et.  al.  (1983)  and  Nair  (1986)  made  pioneering  contribu¬ 
tions  to  this  area  of  research.  Ghosh  (1986)  used  the  search  linear 
models  (see  Srivastava  1975)  to  explain  dispersion  effects  in  factor 
screening  experiments. 

We  first  assume  that  for  the  fitted  model  to  the  data  there  is  no 
significant  lack  of  fit.  We  then  propose  three  sets  of  measures  of  dis¬ 
persion  effects  of  m  factors.  All  three  of  them  are  relevant  in  repli¬ 
cated  factorial  experiments  and  two  of  them  are  applicable  to  unrepli¬ 
cated  factorial  experiments.  The  dispersion  effect  of  a  factor  depends 
on  the  dispersion  at  level  l  and  the  dispersion  at  level  0  of  the 
factor.  The  dispersions  at  levels  1  and  0  based  on  the  least  squares 
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residuals  are  correlated  in  most  situations.  We  introduce  a  method  of 


adjusting  residuals  and  then  calculate  the  dispersions  at  levels  1  and  0 
based  on  these  adjusted  residuals.  The  dispersion  based  on  adjusted 
residuals  at  a  particular  level  of  the  factor  is  uncorrelated  with  the 
dispersions  based  on  residuals  and  adjusted  residuals  at  the  other  level 
of  the  factor. 

The  use  of  the  proposed  measures  of  dispersion  effects  in  a 
factorial  experiment  will  give  a  combination  of  levels  of  control  factors 
which  is  optimum  in  view  of  reducing  the  process  variability  due  to  noise 
in  the  experiment.  The  reduction  of  the  process  variability  due  to  noise 
factors  is  an  important  aspect  in  quality  control  studies. 

2.  Dispersion  Effects 

We  consider  a  2m  factorial  experiment  under  a  completely  randomized 
design.  Let  T(nxn)  be  the  design.  The  columns  of  T  denote  factors  that 
are  controllable  at  their  lower  and  upper  levels.  The  rows  of  T  denote 
runs  or  treatment  combinations.  Runs  are  level  combinations  of  control 
factors  that  are  actually  used  in  the  experiment  to  collect  observations 
or  data.  The  design  T  is  called  an  inner  array  for  m  control  factors. 

The  inner  array  T  chosen  for  the  experiment  will  play  an  important  role 
in  subsequent  discussions.  An  important  objective  in  quality  control 
studies  is  to  evaluate  the  sensitivity  of  the  manufacturing  process  to 
noise.  A  list  of  noise  factors  likely  to  affect  the  process  is  first 
made.  Various  level  combinations  of  noise  factors  that  provide  a  good 
representation  of  noise  are  then  considered.  The  matrix  representation 
of  the  level  combinations  of  noise  factors  is  called  an  outer  array  (see 
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Taguchi  and  Wu  1985).  Suppose  that  there  are  r  (_>  1)  level  combinations 
of  noise  factors  in  the  experiment.  For  every  run  in  the  inner  array  T, 
we  collect  r  observations  corresponding  to  r  level  combinations  of  the 
outer  array.  The  r  observations  for  a  run  in  T  are  called  r  replicated 
observations.  The  variability  in  r  replicated  observations  for  a  run  in 
T  is  attributable  to  process  variability  due  to  noise  in  the  experiment. 
The  case  r  =  1  is  called  the  unreplicated  experiment  and  the  case  r  >  1 
is  called  the  replicated  experiment.  Again,  for  simplicity  equal  repli¬ 
cation  is  considered  for  the  replicated  experiment  and  the  idea  is  easily 
extendable  to  unequal  replications. 

Let  y . .  be  the  j th  observation  for  the  ith  run,  y.  be  the  mean  of 
ij  i 

all  observations  for  the  run  i,  i=l,...,n  and  j=l,...r,  and  N  (=  nr) 
be  the  total  number  of  observations.  The  standard  linear  model  for  the 
experiment  is 

E(y)  -  XB,  (1) 

V(£)  =  o2I,  (2) 

Rank  X  =  p,  (3) 

where  ^(Nxl)  is  the  vector  of  observations  and  y_  =  (y  ^  , . . .  ,y  ^  ; . . .  ; 
ynj . .  )  ’  ,  jKp*l)  is  the  vector  of  factorial  effects  considered  in 

the  experiment,  X(Nxp)  is  a  known  matrix  that  depends  on  the  inner  array 

2  —  1 
T  and  a  is  an  unknown  constant.  We  denote  H  =*  X(X'X)  X'  and  R  = 

(I-H).  The  vectors  jr  =  and  y_  -  y  =  Ry  are  the  vector  of  least  squares 

fitted  values  and  the  vector  of  residuals,  respectively.  The  fitted 

values  for  all  observations  corresponding  to  the  ith  run  are  identical 

and  is  denoted  by  y^,  i  =  l . .  Suppose  that  for  the  fitted  model  to 


^niEWTUtW^: V,7V^TTJV'_'VVW*.rvr1r- 


Che  daCa  there  is  no  significant  lack  of  fit.  The  sum  of  squares  of 
nr  »  y 

error  is  SSE  =  Z  Z  (y..-y.)  »  the  mean  square  of  error  is  MSE  = 
i-1  j  =  l  1J  1 

nr 

(SSE/ (N-p) ) ,  the  sum  of  squares  of  pure  error  is  SSPE  =  Z  Z  (y.  -y ,  ) 

i  =  l  j=l  J 

and  the  mean  square  of  pure  error  is  MSPE  =  (SSPE/n( r-1 ) ) .  Note  that 

2 

both  MSE  and  MSPE  are  measures  of  error  variance  a  .  We  now  take  MSE  and 
MSPE  as  descriptive  measures  of  noise.  We  then  express  MSPE  as  the 
weighted  average  of  (MSPE) ^  and  (MSPE)q,  where  (MSPE)^  is  called  the  con¬ 
tribution  of  the  level  u  (u  =  0,1)  of  the  factor  to  MSPE.  Formal  expres¬ 
sions  of  (MSPE)u,  u  =  0,1  are  given  in  the  next  section.  We  do  the  same 
for  MSE.  Different  levels  of  a  factor  may  contribute  differently  to  MSE 
and  MSPE.  In  general  the  contributions  of  levels  of  a  factor  to  noise 
(measured  by  MSPE  or  MSE)  are  called  the  dispersions  at  levels  of  the 
factor.  The  dispersion  effect  of  a  factor  is  the  ratio  of  the  dispersion 


at  level  1  and  the  dispersion  at  level  0  of  the  factor  (see  Box  and  Meyer 
1986).  If  the  dispersion  effect  is  greater  than  1  then  the  dispersion  at 
level  1  is  more  than  the  dispersion  at  level  0.  We  then  prefer  level  0 
over  level  1  in  terms  of  smaller  dispersion.  Similarly  if  the  dispersion 
effect  is  less  than  1  we  prefer  level  1  over  level  0.  Although  we  use 
the  definition  of  dispersion  effect  presented  in  Box  and  Meyer  (1986), 
one  may  also  take  the  logarithm  of  the  proposed  ratio  or  the  difference 
between  dispersions  at  level  1  and  level  0  of  the  factor.  The  use  of  any 
of  these  definitions  will  give  the  same  conclusion  since  each  of  them 
compares  the  dispersions  at  level  l  and  level  0.  The  main  theme  of  this 
paper  is  to  investigate  the  possible  ways  of  measuring  dispersion  effects 
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of  all  factors.  We  would  like  to  make  it  clear  that  the  proposed 
measures  in  this  paper  are  all  descriptive. 

3.  Measuring  Dispersion  Effects 

We  take  a  single  factor  out  of  m  factors  and  develop  the  methods  of 
measuring  the  dispersion  effect  of  the  chosen  factor.  For  simplicity  of 
the  presentation,  we  do  not  introduce  any  notation  for  the  chosen 
factor.  The  chosen  factor  appears  at  levels  1  and  0  in  the  n  runs  of  the 
inner  array  T.  We  now  introduce  an  indicator  variable  to  identify  the 
runs  at  which  the  factor  appears  at  level  1  and  level  0,  respectively. 
Distinguishing  between  the  level  1  and  the  level  0  runs  for  the  factor 
will  enable  us  to  measure  dispersions  at  level  1  and  level  0  of  the 
factor.  We  define  for  i=l,...n, 

f  l  if  the  level  of  the  factor  in  the  ith  run  is  1, 


0  if  the  level  of  the  factor  in  the  ith  run  is  0. 
3.1.  First  Measure 


We  have 


SSPE 


E  1  'Vyii_yi^  +  1  1  C  1_6i  )  ^  yi  i  "yi 

i  =  l  j=l  1  1  i  =  1  j=l  1  1 


-  1 2 


Not  i 


ice  that  Z  (y..-y.)  is  the  total  (corrected)  sum  of  squares  of  r 
j=l  1J  1 

observations  for  the  ith  run  with  the  degrees  of  freedom  (r-1).  The 

first  component  in  SSPE  corresponds  to  level  I  of  the  factor  and  has 

n 

degrees  of  freedom  f  Z  6,)fr-l).  The  second  component  corresponds  to 

i=  1  1 

n 

level  0  of  the  factor  and  has  degrees  of  freedom  (  Z  ( 1-6  ))fr-l).  The 

i-1  1 

set  of  measures  of  dispersions  at  levels  1  and  0  of  the  factor  are 


J 


S^( 1)  = 


z  E  6Jyij-yi)‘ 

i=l  1=1  J 


(  E  «,)  (r-1 ) 
i  =  l 


.£.  £  (i-^Ky^-yi)2 

s2od)  =  iiUii - : _ ,  (4) 

(  E  (l"«5,))  (r-1) 

1=1  1 

respectively.  The  first  measure  of  the  dispersion  effect  of  the  chosen 

2  2 

factor  is  therefore  S^(1)/Sq(1).  Notice  that  this  measure  is  possible 
only  for  a  replicated  experiment. 

3.2.  Second  Measure 

SSE  =  }.  }.  5i(yij-yJ2  +  £  (i-'Si  )(ytj-yi  )2. 

1=1  j=l  J  i= 1  j=l  J 

Note  that  (y..-y.j  is  the  residual  for  the  jth  observation  on  the  ith  run 
v  ij  l 


and  E  (y, .-y. )  is  the  residual  sum  of  squares  for  r  observations  on  the 

j-1  J 

ith  run.  The  first  component  in  SSE  corresponds  to  level  1  of  the  factor 


and  has  V\  degrees  of  freedom  (d.f.).  The  second  component  corresponds 
to  level  0  of  the  factor  and  has  Vo  degrees  of  freedom  (d.f.)  The  set  of 
measures  of  dispersions  at  levels  1  and  0  of  the  factor  are 


Sj(2)  = 


Z  Z  6  .  f  y . . -y . 

•  i  •  i  i  '  ij  i 
i  =  l  j  =  1  J 


E  E  {  1-6  )  f  y  - 


<2/TS  _  1  =  1  j-1 


1  Myij“yi 
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i 

l 

respectively.  The  second  measure  of  the  dispersion  effect  of  the  chosen 
factor  is  S*(2)/S*(2). 

We  now  denote  y^  as  the  vector  of  observations  corresponding  to  the 
runs  with  the  chosen  factor  at  level  u,  u  3  1,0.  Note  that  consists  of 
y^  and  y^.  Let  be  the  submatrix  of  X  corresponding  to  y^,  ^  be  the 

least  squares  fitted  values  for  y^.  The  elements  ^  and  y^-y^,  u  =  1,0, 

are  linear  functions  of  the  elements  in  y.  It  can  be  seen  that  y  = 

—  -i-U 

X  (X’X)  ^X'y,  u  =  1,0.  We  denote  y  -y  =  r  y,  u  =  1,0.  Then  V  =  Rank, 
u  —  hi  mi  u*-  u 

r  and  V  S 2 ( 2 )  =  y'r’r  y,  u=l,0.  Thus  V  S"(2)  is  the  sum  of  squares  of 
u  u  u  —  u  ir  uu 

the  elements  in  r  y,  u  =  1,0. 

u — 

The  second  measure  of  the  dispersion  effect  was  in  fact  proposed  in 
Box  and  Meyer  (1986)  for  unreplicated  experiments.  We  observe  that  Sj(2) 
and  Sq(2)  are  correlated  under  the  model  ( 1  —  3 ) .  In  the  following  sub¬ 
section  we  present  a  method  of  adjusting  the  above  dispersion  measures 
to  make  them  uncorrelated. 

3.3.  Adjusted  Residuals 

Two  vectors  of  residuals  r  and  r^  at  levels  1  and  0  of  the  factor 
are  generally  correlated  under  the  model  f 1-3).  We  now  present  a  vector 
of  "adjusted  residuals"  at  level  0  of  the  factor,  adjusted  w.r.t  r j _y  so 
that  it  is  uncorrelated  with  r^.  Let  be  a  submatrix  of  r^  so 

that  Rank  rj^=v^*  We  write 


r  =  r  f I-r '  r  r*  ^  r  ■  =  r  —  r  r 1  r  r?  \  r 
0a  0  11  11  llj  11'  0  0  11  nil1  1 


-1 


(6) 


It  can  be  seen  that  Cov  r  y,  r  y 

1 1~  0a— 


)  =  0  and  hence  Cov.rv,  r^  y1  =  0. 

1 —  Oa — 


In  other  words,  r ^y_  and  r^a_y  are  uncorrelated.  We  call  fQa_y  the  vector 
of  "adjusted  residuals"  at  level  0  of  the  factor,  adjusted  w.r.t.  the 
residuals  at  level  1  of  the  factor.  It  can  be  checked  that  Rank  r = 
((N-p)-Vj)  =  V  (say).  Let  r^  be  a  (Vq*n)  submatrix  of  r^  with  Rank 


r01  =  ^0*  wr*-t:e 


r  la  rL^1-r01^r01r01  ^  roJ  r  1 _r  L  r0 1  ^  r0 1  r0 1  r01*  f-7) 


.  i-i. 


Again,  Cov(  r^_,  r  =  0.  In  other  words  r^  and  r^a_y  are  uncorrelated. 


We  call  r  y  the  vector  of  "adjusted  residuals"  at  level  1  of  the  factor, 


adjusted  w.r.t.  the  residuals  at  level  0  of  the  factor.  We  have  Rank 


la 


((N-p)-Vn)  =  Vla(say).  We  have 


r0a^  ~  r0'I_rll^rLlril^  rll'!rC^0 


(-1 


rla^-  =  rl'I_r01^r01r01^  r01^ri^r 


(8) 


The  proof  is  given  in  Theorem  7  in  the  Appendix.  Thus  for  u  =  1,0,  r  _y 
depends  on  y_  only  through  y^  and,  moreover,  Cov  ( r  ^  v., r  ja_y)  =  t.e., 

they  are  uncorrelated  under  f 1 —3 ) .  We  now  present  an  illustrative 


example . 

Example  1 

We  consider  the  example  from  Box  and  Meyer  (1986),  page  20,  and 
Taguchi  and  Wu  (1985),  page  68.  Daniel's  normal  probability  plot  in¬ 
dicates  that,  over  the  ranges  studied,  only  factors  B  and  C  affect 
tensile  location  by  amounts  not  readily  attributed  to  noise  (see  Box  and 
Meyer  1986).  We  now  fit  the  following  standard  linear  model  to  the  data 
Ef  y  x  j  ,x  2 ) )  =  U  +  OjB  +  C,  where  =  0,1,  a.  =  (2x.-l),  u  is  the 


general  mean,  B  and  C  are  the  main  effects  of  the  factors.  We  can  write 

the  above  model  in  the  form  (l-3).  Notice  that  N  =  16,  n  =  4,  p  =  3, 
n  n 

(  Z  <5 . )  =  (  Z  (l-S.))  =  2  for  both  factors  B  and  C.  The  F  value  for  the 
i=l  1  i-1  1 

lack  of  fit  test  under  the  assumption  of  normality  is  .1971  (<  l)  and  we 
therefore  conclude  that  there  is  no  significant  lack  of  fit.  The  inner 


array  T  is  given  by 


-  I"1  1 

L1  0 


0  0 
0  1 


Let  us  now  choose  the  factor  B.  Recall  that  the  vector  of 
observations  y  consists  of  level  1  observations  yj  and  level  0 
observations  yg  of  the  factor  B.  We  have  in  y^  4  observations  on  each  of 
the  runs  (l,l)  and  (l,0)  and  thus  y^  =  (42.4,  42.4,  42.4,  42.5,  44.7, 
45.9,  45.5,  46.5).  We  also  have  in  y_Q  4  observations  on  each  of  runs 
(0,0)  and  (0,1)  and  thus  =  (43.7,  42.2,  43.6,  44.0,  40.2,  40.6,  40.6, 


40.2).  We  now  present  the  matrices 


1  1 


1  1  1 

l  l  1 

l  1  1 

1  1  -1 

1  1  -1 

l  1  -l 

l  l  -1 


1  -1  -1 

1  -1  -1 

1  -1  -l 

1  -1  -1 

1  -1  1 

1  -1  l 

1  -1  1 

1  -1  1 


We  write  y'  =  [y_[:  y^],  X'  =  fxj:X^|,  L{  =  Xp'XpX'y,  yQ  = 

“I  *  * 

Xq(x’X)  X'y^  r  L_y  =  anc*  r qZ  =  Xo'^o’  Can  13,2  checked  that 

X'X  =  161^  and  the  matrices  r^fs^lh)  and  r^Sxlb)  are  given  by 


We  now  obtain 


_J 

16 


rl  = 


i  -H 

16 


(-  -) 

16 

I  -  — 

16 


16 


-)  :  4 


-) 
16  ' 

J 

16 


J 

r.  Ji 

:  i 

3J 

(- 

J1 

L6 

16 

Tb 

16 

: 

j 

:  f- 

-)  : 

I 

3J 

16 

~16 

16 

Tb 

where  J  is  a  (4x4)  matrix  with  all  its  elements  unity.  It  can  be  checked 
that  =  Rank,  r^  =  7  and  Vg  =  Rank  Tq  =  7.  Thus  Vga  =  V^a  =  6.  The 

matrix  r  is  obtained  from  r^  by  deleting  the  last  row.  The  calcula¬ 
tions  of  rga  from  (6)  and  the  vector  of  adjusted  residuals  r^jy  are 
straightforward.  The  calculations  of  r^a  and  the  vector  of  adjusted 

residuals  r,  y  are  similar, 
la— 

3 . 4  Third  Set  of  Measures 

Let  rual  be  a  (v^xn)  submatrix  of  rQa  with  rank  rual  =  Vua , 
u  =  1,0.  We  now  have  the  sura  of  squares  of  the  sets  of  linear 

functions  r  ,y  and  r  ,y  [see  Scheffe  19591  as 
ui—  uai—  1 


ss^uly)  =  y'ruifruiruJ~lruiy  * 

SS(rualy)  =  y'r’al[rualr’al]‘lrualy  , 


with  d.f.  V  and  V  ,  respectively  (u  =  1,0). 


i 

Run 

Z  (yn  -  y,  )2/(r-i) 
j=i  J 

i 

11 

42.425 

42.4875 

.0025  ‘ 

2 

10 

45.650 

45.5875 

.  5700 

3 

00 

43.375 

43.4375 

.6425 

4 

01 

40.400 

40.3  375 

.0533 

Table  2 

Numerical  Values  of  Measures  Of  Dispersion  For  Factors  B  and  C 


Factor 

wmtm 

RB 

m 

ran 

so(3) 

B 

.2863 

.3479 

.2498 

.3027 

.2543 

.3071 

C 

.0279 

.6063 

.0234 

.  5241 

.0329 

.  5286 

Table  3 

Numerical  Values  of  Measures  Of  Dispersion  Effects  For  Factors  B  and  C 


Factor 

Sj(l) 

S  j(  2) 

Sj(  3) 

sj(3) 

Sla(3) 

IBM 

mm 

S0(3) 

Sn  (3^ 

0a 

KH 

S0a(3) 

B 

.8229 

.8252 

.8281 

.7310 

.9323 

.8229 

C 

.0460 

.0542 

.0622 

.0543 

.0528 

.0460 
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4.  Interpretation  and  Application  of  the  Measures 


We  now  discuss  the  rationale  for  the  proposed  measures  ot  dispersion 
effects  and  how  they  can  be  used  in  determining  the  optimum  level  of  each 
control  factor  in  view  of  reducing  the  process  variability.  We  have 


MSPE  = 


/.V1-*  iV 

S?(1)  +  — - n -  SQ 


SqCO 


Thus  S^(l)  and  Sq(1)  are  regarded  as  fMSPE)i  and  (MSPE)o  in  the  notation 

2  2 

of  Section  2.  If  the  first  measure  of  the  dispersion  effect  S^( 1)/Sq( 1) 
is  greater  than  1.  We  then  say  that  level  0  of  the  factor  has  less 

contribution  to  MSPE  and  therefore  would  be  preferred  to  level  1  in 

2  2 

view  of  stability  against  noise  factors.  If  S|( 1)/Sq( 1)  is  less  than 
l,  level  1  would  be  preferred  to  level  0.  We  observe 

MSE  =  ( _ L\  S2  ( 2)  +  ( - ^S2(2). 

y<N-p )j  1  \(N-P )j  0 

If  the  second  measure  of  the  dispersion  effect  S^(2)/S^(2)  is  larger  than 
1,  we  then  conclude  that  level  0  of  the  factor  has  less  contribution  to 

MSE  and  therefore  would  be  preferred  to  level  1  in  view  of  reducing 

2  2 

process  variability  due  to  noise  factors.  If  Sj(2)/S^(2)  is  smaller  than 
1,  level  1  would  be  preferred  to  level  0. 

We  notice  that 

MSE  *  (— L\  S2(3)  +  ( S2  (3)l  =  \( ^L.\  S2  (3)  +  (— 2— \  S2(  3)  . 


The  third  set  of  measures  of  the  dispersion  effect  consists  of  4  measures 
of  the  dispersion  effect.  As  before,  the  numerical  value  of  a  measure 
greater  than  1  indicates  a  preference  of  level  0  over  level  1  and  the 
numerical  value  less  than  1  indicates  a  preference  of  level  1  over  level 
0  in  terms  of  smaller  dispersion.  If  the  numerical  values  of  all  pro¬ 
posed  measures  are  near  1,  then  both  levels  1  and  0  are  equally  prefer¬ 
able  in  terras  of  dispersion. 

Example  1  (continued) 

We  find  from  Table  3  that  numerical  values  of  all  measures  of  dis¬ 
persion  effects  are  less  than  1  for  both  factors  B  and  C.  We  thus  prefer 
level  1  over  level  0  for  both  factors.  The  best  choice  for  level  combi¬ 
nation  of  control  factors  B  and  C  is  therefore  (1,1)  in  view  of  stability 
against  noise.  Notice  that  the  numerical  values  of  all  measures  are  not 
only  greater  but  also  closer  to  1  for  the  factor  B  than  .the  factor  C. 

The  next  best  choice  for  level  combination  of  B  and  C  is  therefore  (0,1) 
in  view  of  stability  against  noise.  Our  choices  for  the  best  and  the 


next  best  level  combinations  are  also  supported  by  the  numerical  values 
r 

r  'i  2 


of  Z  (y..  -  y.l  /fr-l)  in  Table  1.  This  is  of  course  very  natural 

j=r  lJ  1 


because  the  inner  array  T  in  this  example  consists  of  all  runs  of  a  2 
factorial  experiment.  In  practice  we  would  have  a  fractional  factorial 


instead  of  a  complete  factorial  as  the  inner  array  in  most  situations. 

r 

We  would  then  not  be  able  to  calculate  Z  fy..  -  y.)^/;"r-l]  for  all  runs 


to  calculate  Z  <y..  -  y.)/" r-l]  for  al 
j-1  lJ  1 


but  only  for  those  runs  in  the  fractional  factorial  inner  array.  The 
best  choice  of  level  combination  using  the  methods  described  in  this 


paper  may  or  may  not  be  present  in  the  inner  array. 
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5.  Inner  Array  Influence 

We  now  discuss  Che  influence  of  the  inner  array  on  Che  proposed 
measures  of  dispersion  effects.  We  first  consider  the  situation  where 
the  number  of  runs  in  the  inner  array  equals  the  number  of  parameters  in 
0_  of  (1),  or,  in  other  words  n  =  p.  The  inner  array  is  then  a  saturated 
design.  The  class  of  designs  with  n  =  p  includes  the  known  Plackett  and 

Burman  designs  (see  Plackett  and  Burraan  1947).  It  can  be  seen  that  for 

2  2 

an  inner  array  with  n  =  p,  we  have  S^(l)  =  S^(2),  u  *  1,0.  A  general 

2  2 

characterization  of  an  inner  array  for  which  S  (1)  *  S  (2),  u  *  1,0,  and 

u  u 

n  may  or  may  not  be  equal  to  p,  is  available  in  Theorem  1  in  Appendix. 

We  observe  that  for  an  inner  array  with  n  =  p,  the  first  two  measures  of 
dispersion  effects  are  identical. 

We  next  study  the  measures  in  two  extreme  situations:  (1)  Two  vec- 
A  A  » 

tors  of  residuals  y_\  -  yj  -  r  i _y  and  yo  -  yo  “  r0£  at  levels  1  and  0  of 
the  factor  are  uncorrelated,  i.e.,  r^rg  =  0,  (2)  r^  and  v are  com¬ 
pletely  correlated,  i.e.,  rgi  =  Ar i [  for  some  matrix  A.  In  situation 
(1),  we  have  S^( 3)  =  S^(2)  =  S^a(3),  u  =  1,0.  This  in  turn  implies  that 
all  measures  of  dispersion  effects  in  the  third  set  are  identical  and 
identical  to  the  second  measure.  The  proof  is  available  in  Theorem  3  in 
the  Appendix.  We  thus  see  that  in  situation  1  there  is  no  need  for  the 
adjustment  of  residuals  nor  for  the  third  set  of  measures  of  dispersion 
effects.  In  situation  (2),  we  have  r(ga  =  0,  Vga  =  0  and  SS  :  r^  =  0. 
The  proof  is  given  in  Theorem  4  in  the  Appendix.  We  thus  notice  in 
situation  (2)  that  level  1  of  the  factor  makes  all  contribution  to  SSF. 


and  level  0  does  not  make  any  additional  contribution  to  SSK. 


in  case 


VQa  =  V^a  =  0,  we  have  Vq  =  =  (N-p)»  rQl  =  A  11  an<*  A  aons ingular . 


This  is  a  situation  where  levels  0  and  1  have  equal  dispersions.  This  is 


purely  due  to  the  inner  array  influence.  We  have  V  >  f  E  6.  ;(r-l)  and 

A  <1  .  «  V 

l 


n  n 

V  >  (  E  ( 1-6  ) ) ( r  —  1 )  for  all  inner  arrays.  If  V  =  (  E  6.  ) ( r  —  1 )  then 
Ua  i-1  1  13  i-1  1 


S  i  a  ( 3 )  =  SjCl)  and  VQa  =  (  E  (l-S^Ur-l)  then  SQa(3)  =  The  ptoof 


is  given  in  Theorem  2  in  Appendix.  We  note  that  and  V,  are  both 


nonzero  for  r  >  1.  (We  assume  naturally  that  there  is  at  least 


6 ^  =  1  and  at  least  one  ( 1-6^ )  *  1.)  For  the  case  r  *  1,  at  least  one  of 


V.  and  V,  could  be  zero  or  both  of  then  could  be  nonzero. 
0a  la 


We  now  present  an  example  of  an  inner  array  for  which  S;( 2)  and 


Sq(2)  are  uncorrelated  for  one  factor  but  Sj (2)  and  S^(2)  are  correlated 


for  all  other  factors.  This  example  is  remarkable  in  displaying 


contrasting  influence  of  the  inner  array  in  measuring  S^(2),  u  =  1,0,  and 


in  the  need  for  adjustment  of  residuals  for  different  factors. 


Example  2. 


We  consider  a  2  factorial  experiment,  i.e.,  m  =  5.  We  thus  have  5 


controlled  factors  each  at  2  levels.  Let  the  inner  array  T  (8*5)  be 


22 


We  consider  a  situation  where  the  model  for  a  main  effect  plan  fits  the 


data  adequately,  i.e.,  there  is  no  significant  lack  of  fit.  Therefore 
n  ■  8  and  p  =  6.  The  first  column  of  X  has  all  entries  unity.  The 
distinct  rows  of  the  remaining  columns  of  X  are  obtained  from  T  replacing 
0  by  (-1)  and  each  distinct  row  is  replicated  r  times.  We  denote  an 
(r*r)  matrix  with  all  elements  unity  by  J  and 


It  can  be  easily  seen  that 


It  now  follows  that  rjtg  =  0  for  the  factor  3  but  tjtg  *  0  for  factors 
1,2,4  and  5.  Thus  S^(2)  and  ( 2 )  are  uncorrelated  for  the  factor  3  and 
are  correlated  for  factors  1,2,4  and  5. 
b.  Properties 

We  now  state  some  properties  of  the  descriptive  measures  under  the 

2 

model  (1-3).  We  first  observe  that  the  measures  of  dispersion  S  ^ (  1 )  and 
Sq(1)  do  not  depend  on  the  fitted  model  and  all  other  measures  of  dis¬ 
persion  depend  on  the  fitted  model.  The  measures  S^(l)  and  S^(l)  are 

/  2 

always  uncorrelated  under  the  model  (1-3).  The  measures  S  (2)  and  S  (2) 
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nay  however  be  correlated.  They  are  uncorrelated  If  and  only  if 

2  2 

r^r^  =  0.  Similarly  S^(3)  and  Sg(3)  are  generally  correlated  and  they 

are  uncorrelated  if  and  only  if  r^r^  =  0,  or,  equivalently  r^r^  =  0. 

The  dispersion  measures  S2(3)  and  S2  (3),  S2  and  S2(3),  S2  (3)  and 

i  Oa  la  U  la 

2 

Sq  (3)  are  all  uncorrelated  under  (1-3).  Under  the  assumption  that 
£  -  N(X0,  o2l) ,  the  measures  of  the  dispersion  effect  S2(I)/Sq(1), 

S2(3)/S2  (3),  S2  (3)/S2(3)  and  S2  (3)/S2  (3)  have  the  central  F  distri- 

1  0a  la  0  la  0a 

2  2 

bution  with  appropriate  degrees  of  freedom.  The  measures  S  (2)/S  (2)  and 

2  2 

S  (3)/S  (3)  have  the  central  F  distribution  if  and  only  if  r^r^  *  0.  We 

2  2 

question  the  use  of  measures  of  the  dispersion  effect  S^(2)/Sq  (2)  and 

2  2 

S1(3)/Sq(3)  unless  r^r^  =  0.  We  of  course  realize  that  the  condition 
r^rg  =  0  is  too  stringent  to  satisfy  even  for  one  out  of  m  factors. 

We  know  that  the  sum  of  least  squares  residuals  is  equal  to  zero. 

We  observe  here  that  the  suras  of  the  residuals  at  levels  1  and  0  for  the 
chosen  factor  are  zero.  Moreover  the  suras  of  the  adjusted  residuals  at 
levels  1  (adjusted  for  level  0)  and  0  (adjusted  for  level  1)  are  also 
zero.  The  proof  is  given  in  Theorem  6  in  Appendix.  These  observations 
are  useful  in  calculating  the  degrees  of  freedom. 


£ 


7.  Conclusions 


In  Industrial  experiments  for  quality  improvement,  dispersion 
effects  of  factors  play  an  important  role.  They  are  instrumental  in  the 
choice  of  an  optimum  combination  of  levels  of  control  factors.  This 
article  presents  the  descriptive  methods  of  measuring  dispersions  and 
dispersion  effects  at  the  preliminary  stage  of  investigation.  The 
outcome  of  such  comparisons  will  suggest  more  appropriate  complex  models 
for  further  investigation.  We  however  believe  that  the  implementation  of 
the  methods  discussed  in  this  article  will  result  in  highly  informative 
conclusions.  Although  in  this  paper  we  find  the  measures  of  dispersion 
at  levels  of  a  factor  using  one  factor  at  a  time,  the  same  approach  can 
be  used  to  find  the  measures  of  dispersion  at  level  combinations  of 
factors  using  many  factors  at  a  time.  Unless  the  number  of  observations 
is  sufficiently  large  in  every  cell,  the  reliability  of  the  measures 
using  many  factors  will  be  questionable. 
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APPENDIX 


We  now  present  the  proofs  of  many  statements  that  we  have  made  in 


the  main  body  of  the  paper.  We  also  present  some  valuable  technical 


results  in  the  investigation. 


Let  Dj(NxN)  be  a  diagonal  matrix  with  n  sets  of  diagonal  elements 


and  the  elements  in  the  ith  (i  =  l,...n)  set  are  equal  to  6^.  We  define 


Dg  *  I  -  Dj.  It  can  be  seen  that  DiDg  =  0  and  both  Dj  and  Dg  are 


iderapotent  matrices.  We  have  R  =  D^R  +  DgR.  The  matrices  r^  and  rg 


defined  in  Section  3  are  in  fact  non-null  row  vectors  of  DjR  and  DgR, 


respectively.  It  can  be  seen  that  RD^R  *  r^ru»  vu  S^(2)  =  ^'RD^Rv,  and 


sse  =  =  y ,RDiRy  +  y ,RDgRx.  u  =  i»o. 


2  2 

We  now  investigate  the  situation  where  S  (1)  =  S  (2),  u  =  1,0.  In 

u  u 


2  2 

other  words,  we  characterize  the  inner  arrays  for  which  S  (1)  =  S  (2), 

u  u 


u  =  1,0.  We  denote  the  row  of  the  matrix  X  corresponding  to  the  run  i  by 


>cl(lxp).  Note  that  for  each  i,  i  =  l,...,n,  the  row  si  is  repeated  r 


times  in  X.  Let  X*(nxp)  be  a  matrix  whose  ith  row  is  _x ! .  Notice  that 


rows  of  X*  are  in  fact  distinct  rows  of  X.  We  have  X'X  =  r(X*'X*). 


Theorem  1.  We  have  S^(l)  =  ( 2 ) ,  u  =  1,0  if  and  only  if  X*(X*'X*)_1 


X*'  =  in. 


Proof.  Note  that  S^(  1)  =  S^(2),  u  =  1,0,  hold  if  and  only  if  y^  =  y^.  , 


i  =  l,...,n.  The  condition  y^  =  y^,  i  =  l,...,n,  holds  if  and  only  if 


-1  _  for  i.  =  i. 

x!  (X'X)  x!  =  <  r 

~U  _12  1 


0  for  ij  *  121  i 1 ,  i 2 E { 1 , • • • ,n  } . 


The  above  condition  may  be  expressed  as  X*(XfX)  X*  =  —  1^,  or, 

equivalently,  X*(X*'X*)~^  X*'  =  I  .  This  completes  the  proof. 

n 

When  n  =  p,  X*(nxn)  satisfies  the  condition  X*(X*’X*)  3  X*'  =  1^. 

Therefore,  for  n  =  p,  we  get  S2(l)  =  S2(2). 

u  u 

We  now  present  results  showing  the  influence  of  the  inner  array  on 

the  measures  of  dispersion  S2(3),  S2  (3),  S 2  ( 2 )  and  S2(l'),  u  =  1,0. 

K  u  ’  ua  u  u 

Theorem  2. 


1.  vla  2  (  T  ^Hr-1),  VQa  >  i  £  (l-6.))(r-l)  , 
i=l  i=l 


ii.  VuS2la(3)  >  (  £  6ikr-l)S;(l), 
i  =  1 


V0aS0a<3>  2  (ifiO-«i))(r-l)S0(l)  . 


i  i  t .  If  V, 


(  £  6.}(r-l)  then  S^a( 3)  =  S2(l)  , 
i  =  l 


iv.  If  VQa  =  (  Z  (l-6.)](r-l)  then  S2g(3)  =  S2(i). 
i  =  1 

Proof.  It  can  be  checked  that  Cov ; y  .  .  ,  y  -y  )  =  Covfy.,  y  -y  )  and 
-  ij  u  u'  '■'i’  yu  Ju; 

therefore  Cov'y..-v.,  y  -y  1  =  0.  Moreover,  Cov I y . . -y . ,  y  -y  I  =  0 , 
ij  '  i  u  u'  ij  i  Juw  yu 

i  *  u.  It  now  follows  that  any  contrast  of  (y. . -y.  1,  i  =  1 . r  for 

ij  i 

a  fixed  i  with  6.  =  1,  is  orthogonal  to  any  contrast  of  C y  -y  1, 
i  uw  u 

w  •=  l,...,r  for  a  fixed  u  with  I  1-6  ]  =  Furthermore,  any  contrast  of 

^  y . — y .  1  for  all  i  with  6.  =  1  is  orthogonal  to  any  contrast  of  fy  -v  1, 
w  =  l,...,r,  for  a  fixed  u  with  (  1-6  )  -  1.  The  results  (i-iv)  follow 


immediately  from  the  above  facts,  the  relationship  between  the  rank  and 


the  number  of  orthogonal  contrasts,  and  the  fact  that  the  sum  of  squares 
is  equal  to  the  sun  of  suns  of  squares  of  orthogonal  contrasts.  We  now 
study  the  measures  in  two  extreme  situations:  (i)  r ^ y_  and  r  ^y_  are  un¬ 
correlated,  i.e.,  r^rg  =  0>  (ii)  r  ^y_  and  r^  are  completely  correlated, 
i.e.,  r^  =  Ar  ^  for  .some  matrix  A. 

Theorem  3.  Consider  the  situation  r  r'  =  0.  Then  S2(3)  =  S 2 ( 2 )  = 
-  10  u  u 

S2  (3),  u  =  1,0. 
ua 

Proof .  We  first  show  that  S 2  ( 3 )  =  S2(2),  or,  in  other  words,  SSfr^yJ  = 
y_'  r  ir  We  observe  that 


rii  [rnrli 


rn  z 


y'  rj 


f  r  r  ’ 

1  11  11 


,-l 


rll 


0 

0  0 
~  —  r[  Z>  where  r^rj  =  r ^ 

=  r!ri  1  =  (rlri  ^rlrl  £ 
=  (r^)’  =  Z.’riril- 


It  follows  from  the  representations  of  MSE  in  terms  of  (s^(2),  Sq(2))  and 


C  3  ^ ( 3 )  and  Sq3 ( 3 ) )  chat  =  VqS^(2).  The  condition  r^r q  =  0 

implies  that  V^a  =  Vg.  Thus  Sga(3)  =  Sq(2).  The  rest  is  similar.  This 
completes  the  proof. 

Theorem  A.  If  r. 


01 


Ar  then  we  have  rga  =  0,  Vga  =  0  and 


SS|,rQal.zJ  °* 


Proof .  We  write  r^  =  A^r ^  for  a  matrix  Ag  whose  independent  rows  are 


rows  of  A.  This  implies  that  rgrj^  =  ^Qriir'ii  an<^  thus 


0.  The  rest  is 


A0  =  r0riprllrll  *  Hence  from  (6)  we  §eC  r0a  = 

clear.  This  completes  the  proof. 

We  now  present  some  results  which  are  useful  in  establishing 
properties  of  the  proposed  measures  in  the  paper. 

Theorem  5.  Suppose  y_  N(xj?,o  i).  A  necessary  and  sufficient  condition 
that 


y'r[riX  2 

(1)  _  ~  central  y  with  d.f. 

o 

2 

(2)  - ^ central  y  with  d.f. 


Trace  rjr^, 


Trace 


r0r0* 


a 

(3)  and  furthermore,  (1)  and  (2)  are  statistically  independent,  is  that 
rlr0  =  °- 

It  can  be  seen  (see  Rao  1973)  that  a  necessary  and  sufficient  condition 

of  (1),  (2)  and  (3)  to  be  true  is  that  rJr^rgrQ  =  0.  The  condition  is 
equivalent  to  r^r^  =  0* 

Theorem  6. 


3  • 


b. 


j_'  =  0  and  j_'rua  =  0  ,  u  =  1,0,  where  j_'  is  a  vector  with  all 

elements  unity, 

If  r^rg  =  0,  then  for  u  =  1,0, 

b.l.  r  r'  is  an  idempotent  matrix, 
u  u 


b .  2. 


y  -  y  , 

,j!-u  ^-u 


=  r  r  y 

U  I*M1 


b.  3.  X'  r  r’  =  0  , 
u  u  u 

n  r  *  . 

b.4.  z  z  5 . y .  f  y . .  -  y . 

.  ,  .  ,  1  1  1J  l 

1=1  J=1 


n 

=  Z 


z  ( i  ~  6, )y. (y. .  -  y. 1  =  o. 

i=l  j=l  J  1 


Proof.  The  result  (a)  follows  by  considering  the  columns  of  X  for  the 
general  mean  and  the  factor  chosen  and  from  the  fact  that  X'R  =  0.  The 
results  b.l  and  b.2  follow  directly  from  the  structure  of  R  and  the  fact 
that  R  is  an  idempotent  matrix.  The  result  b.3  follows  from  X'R  =  0. 
From  b.3,  we  get  B'X'r  r'y  =  0,  i.e.,  y'r  r'y  =  0.  The  result  b.2 

implies  the  ^(y^  ~  )  =  0  and  hence  the  result  b.4  is  true.  This 

completes  the  proof. 

Theorem  7.  For  Tq3  and  r^  in  (6)  and  (7),  the  equation  (8)  holds. 

Proof .  We  prove  the  equation  (8)  for  rQay.'  Pro°f  for  '•s 

similar.  The  fact  RX  =  0  implies  that  r^y  =  0,  r^z  =  0  ai1d  hence 

T Qay_  ~  0*  Since  z  “  X  =  ^X  an<^  =  we  have  X  =  X  +  R*Z  =  Z  +  rlZ[  + 
r^y^.  The  independent  rows  of  r  ^  are  in  r^  and  the  other  rows  of  r  ^ 
can  be  written  as  Qr ^  for  some  matrix  Q.  The  fact 

U  -  r[i^rnriJ  1  r 1 1 1 1  =  0  imPUes  ( 1  -  rii^rLirii)  1  r 1 1 ^ r [  =  0  and 

hence  r^rj  =  0.  Hence  r^y  =  r^r^. 


This  completes  the  proof. 


